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Vibron-polaron in a-helices. II. Two-vibron bound states 
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The two-vibron dynamics associated to amide-I vibrations in a 3D a-helix is described according to 
a generalized Davydov model. The helix is modeled by three spines of hydrogen-bonded peptide units 
linked via covalent bonds. It is shown that the two-vibron energy spectrum supports both a two- 
vibron free states continuum and two kinds of bound states, called TVBS-I and TVBS-II, connected 
to the trapping of two vibrons onto the same amide-I mode and onto two nearest neighbor amide-I 
modes belonging to the same spine, respectively. At low temperature, non vanishing interspine 
hopping constants yield a three dimensional nature of both TVBS-I and TVBS-II which the wave 
functions extend over the three spines of the helix. At biological temperature, the pairs are confined 
in a given spine and exhibit the same features as the bound states described within a one-dimensional 
model. The interplay between the temperature and the 3D nature of the helix is also responsible 
for the occurrence of a third bound state called TVBS-III which refers to the trapping of two 
vibrons onto two different spines. The experimental signature of the existence of bound states is 
discussed through the simulation of their infrared pump-probe spectroscopic response. Finally, the 
fundamental question of the breather-like behavior of two-vibron bound states is addressed. 

PACS numbers: 63.20.Dj,63.22.+m,71.38.Ht,87.10.+c 



I. INTRODUCTION 

In low-dimensional molecular lattices, the nonlinear 
nature of vibrational excitons (vibrons) plays a key role 
for energy transfer as well as energy storage in both phys- 
ical, chemical and biological systems. This feature has 
become more and more central during the last decade 
through the concept of discrete breathers 0,0,11] . In clas- 
sical anharmonic lattices, discrete breathers correspond 
to time-periodic and spatially localized solutions which 
result from the interplay between the discreteness and 
the nonlinearity. These solutions do not require integra- 
bility for their existence and stability and it has been 
suggested that they should correspond to quite general 
and robust solutions. Since discrete breathers yield a lo- 
cal accumulation of the vibrational energy which might 
be pinned in the lattice or may travel through it, they 
are expected to be of fundamental importance. 

At present, because the occurrence of classical 
breathers is a relatively well understood phenomena, a 
great attention is done to characterize their quantum 
equivalent for which less numerous results are known 
j4j. In that context, bound states involving two high fre- 
quency vibrational excitons, also called two-vibron bound 
states (TVBS), play a central role. When two vibrons are 
excited, the vibrational anharmonicity is responsible for 
the trapping of the two quanta over only a few sites with 
a resulting energy which is lesser than the energy of two 
quanta lying far apart. The lateral interaction yields a 
motion of such a state from one site to another, thus 
leading to the occurrence of a delocalized wave packet 
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with a well-defined momentum. As a result, TVBS are 
the first quantum states which experience the nonlinear- 
ity and can thus be viewed as the quantum counterpart 
of breathers or soliton excitations [a, IS II] • 

The present paper is thus devoted to the characteri- 
zation of the TVBS associated to the amide-I vibrations 
in a 3D a-helix. Since the pionneer work of Davydov 
[1, , the delocalization of amide-I vibrations yields the 
occurrence of vibrons which are expected to play a funda- 
mental role for the transport of the energy released by the 
hydrolysis of adenosine triphophate [Til Il2| . More pre- 
cisely, the energy transfer is achieved by small polarons 
which originate in the strong couplin g b etween the vi- 
brons and the phonons of the helix [ mull ull ITfi It7| . A 
small polaron corresponds to a vibron dressed by a vir- 
tual cloud of phonons describing a localized lattice dis- 
tortion which follows instantaneously the vibron. 

In the paper I of this series, the single-vibron states 
in a 3D model of an a-helix |l8j have been character- 
ized. It has been shown that the vibron dynamics re- 
sults from the competition between two kinds of hop- 
ping processes. The first kind refers to vibrational tran- 
sitions between different spines whereas the second kind 
involves vibron hops along the same spine. This study 
has revealed that several parameters such as the tem- 
perature, the small polaron binding energy and the helix 
backbone conformation, allow for a transition between 
two regimes. At low temperature or weak small polaron 
binding energy, the polaron behaves as an undressed vi- 
bron delocalized between the different spines. By con- 
trast, at biological temperature or strong small polaron 
binding energy, the dressing effect strongly reduces the 
vibrational exchanges between different spines so that the 
polaron tends to propagate along a single spine, only. Al- 
though the phonon spectrum exhibits both an acoustic 
branch and an optical branch, it has been shown that the 
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previous features originate in the coupling between the 
vibrons and the acoustic phonons. 

In recent works [lfi l2fT |. TVBS have been studied 
within the ID model of an a-helix in which the vibron 
dynamics is reduced to that of a single spine of hydrogen- 
bound peptide unit. In that case, it has been shown that 
both the intramolecular anharmonicity of each amide-I 
vibration and the strong vibron-phonon interaction act as 
nonlinear sources which break the vibron independence 
and favor the formation of two kinds of bound states. 
In the harmonic situation, the two bound states appear 
as combinations of states involving the trapping of the 
two vibrons onto the same amide-I mode and onto near- 
est neighbor amide-I modes in a given spine. By con- 
trast, the intramolecular anharmonicity reduces the hy- 
bridization between these two kinds of trapping so that 
the low frequency bound state refers to the trapping of 
the two vibrons onto the same amide-I mode whereas 
the high frequency bound state characterizes the trap- 
ping onto nearest neighbor amide-I vibrations in a given 
spine. These results were corroborated by the pump- 
probe spectroscopy experiment of the N-H mode in a 
stable a-helix which has revealed the two excited state 
absorption bands connected to the two kinds of bound 
states |2lll22|. 

The present paper is organized as follows. In Sec. II, 
the model to describe the vibron-phonon dynamics in 3D 
a-helices is introduced. Then, the procedure to remove 
the intramolecular anharmonicity and to renormalize the 
vibron-phonon interaction in order to obtain the dressed 
anharmonic vibron Hamiltonian is summarized. Finally, 
the number states method used to solve the two-vibron 
Schrodinger equation is outlined. In Sec. Ill, a detailed 
analysis of the two-vibron energy spectrum is performed 
depending on the values taken by the relevant parameters 
of the problem. The results are discussed and interpreted 
in Sec. IV. 



II. THEORETICAL BACKGROUND 

As pointed out in the introduction, the two-vibron dy- 
namics is strongly sensitive to the nonlinear nature of 
the vibron-phonon system. In a-helices, this nonlincar- 
ity originates from both the intramolecular anharmonic- 
ity of each amide-I vibration and from the strong vibron- 
phonon coupling . In the present section, we thus briefly 
summarize the model described in paper I |18| and in- 
troduce the anharmonic nature of each amide-I mode ac- 
cording to the ID model detailed in Ref. [l9j. 

A. Helix structure and model Hamiltonian 

Let us consider a sequence of TV amino acid units 
(called residues) regularly distributed along a polypep- 
tide chain and let n the index which labels the position 
of the nth residue. In three dimension, the structure of 



the polypeptide chain is stabilized by the hydrogen bonds 
between the residues n and n + 3. The resulting confor- 
mation is a 3D a-helix in which each residue is related 
to the next one by a translation of h = 1.5 A and by a 
rotation of 6q = 100° leading to 3.6 residues per turn of 
helix. The N residues are assumed to behave as point- 
like entities which the equilibrium positions are located 
on sites distributed along the helix which the radius is 
fixed to Rq — 2.8 A The geometry of the helix is 

thus specified in the (xyz) cartesian frame where the z 
direction is parallel to the axis of the helix (see Fig. la 
in paper I 0). 

The nth site contains a amide-I vibration which be- 
haves as an internal high frequency oscillator described 
by the standard creation and annihilation vibron opera- 
tors and b n . By introducing the intramolecular anhar- 
monicity of each amide-I mode, the vibron Hamiltonian 
is written as (within the convention h=l) 

h v = Y / uj obib n +j 3 (bi + b n ) a +^(bi + b n f + ... 

n 

- 5EE J ( n - n ')(&! + M(^ + M (!) 

n n' 

where luq stands for the internal frequency of the nth 
amide-I mode and where J(n — n') denotes the lateral 
hopping constant between the residues n and nl . In 
Eq.Q, 73 and 74 represent the cubic and quartic an- 
harmonic parameters of each amide-I mode. 

The amide-I vibrations interact with the phonons of 
the helix which corresponds to the dynamics of the ex- 
ternal motions of the residues. As detailed in paper 
I 01, this dynamics is expressed in N local frames 
(e r (n),eg(n),ejn)) attached to each residue (see Fig. 
lb in paper I [lg) so that the displacement of the nth 
residue in the nth local frame corresponds to a Bloch 
wave with wave vector k |24| . Therefore, the normal 
mode decomposition is achieved by performing the di- 
agonalization of a (3 x 3) dynamical matrix D(k) for 
each k value. Such a procedure allows us to define three 
eigenvalues £l\ a and three eigenvectors e^ G labeled by the 
index a — 1,2,3 associated to the 3D nature of the he- 
lix. The two indexes k and a specify a particular phonon 
mode with energy Mlka, quasi- momentum hk and po- 
larization efe CT . The quantum dynamics of each mode is 
described by the creation at and annihilation a ka opera- 
tors so that the phonon Hamiltonian is written as (within 
the convention h = 1) 

H P = J2 n k*(al a a ktr +±) (2) 

ka 

Finally, according to the Davydov model, the vibron- 
phonon interaction originates in the modulation of the 
vibrational frequency of each amide-I vibration by the 
external motion of the residues. Within the deforma- 
tion potential approximation, the vibron-phonon cou- 
pling Hamiltonian is written in terms of the phonon nor- 
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mal mode coordinates as 



correction n = 120(73/wo) 2 — 1274/wo as 
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where A kai which accounts for the modulation of the 
frequency of the nth amide-I vibration due to its coupling 
with the phonon mode specified by the wave vector k and 
the index a, is expressed as 



A ka = 



X\m\C k (m)e* k0 



(4) 



In Eq.Q, X\n-n'\ denotes the strength of the interaction 
between the nth amide-I mode and the external displace- 
ment of the n'th residue and Cfc(m) is a vector which 
account for the helix geometry and which is defined in 
Eq. (16) of paper I [3. 



B. Effective Hamiltonian for the two-vibron 
dynamics 

The vibron-phonon dynamics is described by the full 
H„ + AH vp which cannot be 



Hamiltonian H = H v 



solved exactly due to the nonlinear sources. Neverthe- 
less, as detailed in Ref. |l9j and outlined in paper I, it 
can be expressed in an improved way by performing two 
unitary transformations. First, by following Kimball et 
al. |fj, a unitary transformation is applied to remove 
the intramolecular anharmonicity of each amide-I mode. 
Then, a modified Lang-Firsov transformation is realized 
to renormalize a part of the vibron-phonon interaction 
[25|. Finally, a mean field procedure is applied to sepa- 
rate the vibron degrees of freedom from the phonon co- 
ordinates in the residual coupling so that the effective 
dressed anharmonic vibron Hamiltonian is written as 

H eff =J2"ob+b n -Ab+ 2 bl 

n 

-J2lz(n-n')b+b+b n b n , 

7in' 

-Y^Mn- n')$ nn ,(N n + N n ,)b+b n , 

11 a' 

- £ Mn - n')$ nn ,(N n + N n ,fbt 2 bl> 

nn' 

- Y Mn - n')$ nn >(N n + N n ,)b+[N n + N n ,]b n ,(5) 



where the different parameters enter Eq.JSJ) are ex- 
pressed in terms of the intramolecular anharmonicity 
A = 3O7 2 /cl>o — 674, the generalized small polaron binding 

IA I 2 

energy £ , s(n) = ^2 ka . ^ cos(fcn) and the anharmonic 



(b = u -2A-Y, n B{n)/2-{l + Arj)E B (0) 

A = A+(l + 8n)E B (Q) 

B{n) = B(n) + 2(1 + 4n)£ s (") 

B{n) = 144 J(n) (73 /w„) 2 

Ji(n) = J(n)(l + 44( 73 /wo) 2 - 1274/wo) 

J 2 (n) = 4J(n)( 73 M ) ) 2 

J 3 (n) = J(n) (22(73 M)) 2 ~ 1274/wo) (6) 

In Eq.©, <& nn >{X) =exp(-S , (n-n / ,T)[l + 2n + 2nX]) 
where S(n, T) is the coupling constant defined as (ks 
denotes the Boltzmann constant) 



S(n,T) = Y 



ka 



n kr7 1 



coth( Sl )(1 " cos(M) (7) 



As clearly shown in Eq.©, the nonlinear sources 
strongly modify the dynamics of the dressed anharmonic 
vibrons. They are first responsible for a redshift of each 
amide-I frequency, i.e. loq — > ujq. Then, they lead to the 
occurrence of coupling terms which break the vibron in- 
dependence. The terms Ab^ 2 b„ and B(n — n')b+b~^,b n b n i 
yield an attractive interaction between two vibrons and 
favor their trapping around neighboring amide-I sites. 
Finally, the nonlinearities modify the vibron hopping pro- 
cesses via the dressing effect which reduces the hopping 
constants connecting two sites n and n' according to the 
dressing function . Note that the contributions pro- 
portional to J 2 characterize hops in the course of which 
two vibrons realize simultaneously a transition whereas 
the terms proportional to J3 affect single-vibron hops 
from states formed by two vibrons located onto the same 
amide-I vibration. These two latter mechanism are spe- 
cific of the two-vibron dynamics and do not occur in the 
single-vibron dynamics. At this step let us mention that 
the 3D nature of the helix provides a site dependent small 
polaron binding energy Esin) as well as a site dependent 
coupling constant S(n — n', T). As it will be discussed in 
the following section, these two features discriminate be- 
tween interspine and intraspine processes for both vibron 
hops and vibron-vibron interactions. 

To characterize the two-vibron dynamics, the corre- 
sponding Schrodinger equation H e ff \ ty) = ui | has 
to be solved. Since the Hamiltonian H e ff (Eq.©) con- 
serves the number of vibrons, this can be achieved by 
using the number states method § which was success- 
fully applied to molecular adsorbates [26Tl2?ll2^| and ID 
model of a-helices |2(j . Within this method, the two- 
vibron wave function is expanded as 



m = 



^ *(ni,n 2 ) I ni,n 2 ) 



(8) 



ni ,n2>ni 



where {| ni,n2)} denotes a local basis set normalized 
and symmetrized according to the restricting n 2 > n\ 
and where a particular vector | ni,n 2 ) characterizes two 
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vibrons located onto the sites n\ and ri2, respectively. 
This basis set generates the entire two-vibron subspace 
which the dimension N(N + l)/2 represents the number 
of ways for distributing two indistinguishable quanta onto 
N sites. 

The Schrodinger equation can be expressed in an im- 
proved way by taking advantage of the helix symmetry. 
Indeed, as for the phonon dynamics, the two-vibron wave 
function is invariant under a rotation by an angle of 6q 
around the helix axis followed by a translation along this 
axis by a distance h. As a result, it can be expanded as 
a Bloch wave as 

tf(ni,na = rn+rn) = ^= ^ e iK <- ni+m / 2 ^ K (m) (9) 

where the total momentum K, which takes N values be- 
longing to the first Brillouin zone of the helix, is asso- 
ciated to the motion of the center of mass of the two 
vibrons. Therefore, the resulting wave function <J if (m) 
refers to the degree of freedom m which characterizes 
the distance between the two vibrons along the helix. 
Since the momentum if is a good quantum number, the 
Hamiltonian H e ff appears as block diagonal and the 
Schrodinger equation can be solved for each K value. 
This procedure, achieved numerically, allows us to de- 
fined a set of two-vibron eigenenergies u>k\ an d a set of 
two-vibron eigenstates ^ k\(iti) which are specified by 
the index A = 1, (N + l)/2 for each K value. 

In the following section, this procedure is illustrated 
to characterize the two-vibron energy spectrum of a 3D 
a-helix. 



III. NUMERICAL RESULTS 

To describe the vibron dynamics in 3D a-helices, the 
hopping constant al ong the hydrogen bonds is fixed to 
J(3) = 7.8 cm -1 [ill H(| whereas the hopping con- 
stants between different spines of hydrogen-bonded pep- 
tide units, calculated in Ref. j2^ are equal to J(l) = 
— 12.4 cm -1 and J(2) = 3.9 cm -1 . The other elements 
of the hopping matrix, rather small, are listed in the Ta- 
ble I of Ref. [2^| . The harmonic frequency is fixed to 
loq = 1695 cm -1 and the anharmonic constant is equal 
to A = 8.0 cm" 1 [H HJ UK. The cubic and quartic 
anharmonic parameters can be expressed approximately 
by using the relation Wj^/uio w 674 w A |l9j. Note that 
these parameter values are in a rather good agreement 
with recent ab initio calculations |3^ |. 

As discussed in paper I E3 , the harmonic dynamics of 
the phonons is essentially governed by the force constant 
Kv n _ n i\ which refers to the second derivative of the pair 
potential connecting the two residues located onto the 
sites n and n' . The force constant K2 is set to zero and 
we assume Ki n i = for n > 4. The force constant of the 
hydrogen bonds is fixed to K 3 = 15 Nm" 1 [H El El US 
l23l whereas the force constant of the covalent bonds is 



equal to K\ = 60 Nm" 1 23]. Finally, the mass M which 
enters the phonon dynamics has been fixed to 2.0 10~ 25 
kg. The strength of the vibron-phonon coupling Eq.lQ 
is given by the parameters Xln-n'l UM which accounts 
for the modulation of the nth amide-I frequency due to 
the external motion of the n'th residue. According to the 
ID Davydov model, the parameter \3 ranges between 35 
and 62 pN and we choose X2 = whereas we treat xi as 
a parameter smaller than \3 |23j. Note that Xn = for 
n > 4. 

The two-vibron energy spectrum is shown in Figs. 1 
for T = 5 K, xi = 25 pN and X 3 = 50 pN (E B (0) = 8.9 
cm -1 ). The spectrum, centered onto the corrected fre- 
quency 2u!q (Eq.@)), describes the two-vibron dispersion 
curves drawn in half of the first Brillouin zone of the he- 
lix, i.e. < K < 7r. In the harmonic approximation, 
i.e. A = (Fig. la), the spectrum exhibits an energy 
continuum located around 2a>o with a bandwidth equal 
to 101.9 cm . This continuum contains the states de- 
scribing two independent vibrons and called two-vibron 
free states (TVFS) El- 

For a given K value, a free state 
refers to two independent vibrons with wave vectors q\ 
and (72 and energy u> qi and u q2 so that K = qi + 92 and 
luk\ = +^<j2- Below the continuum, the spectrum ex- 
hibits two bands connected to two different bound states. 
The low frequency band, located below the TVFS over 
the entire Brillouin zone, refers to bound states called 
TVBS-I. The binding energy of TVBS-I, i.e. the gap be- 
tween the zero wave vector TVBS-I and the bottom of the 
TVFS continuum, is equal to 4.8 cm -1 whereas its band- 
width is about 24.0 cm -1 . The second band, which refers 
to bound states denoted TVBS-II, lies below the con- 
tinuum in a wave vector window which ranges between 
K = 0.5 and K — 0.9. When the anharmonic parameter 
is set to A = 8 cm -1 (Fig. lb), the shape of the TVFS 
continuum does not exhibit any significant modification. 
Note that the center of the continuum is redshifted due 
to the dependence of 2ujq (Eq.©) on the anharmonic 
parameter A (not drawn in the figure). Nevertheless, 
the TVFS bandwidth is slightly reduced to 100.0 cm -1 
so that the wave vector window for the occurrence of 
the TVBS-II opens and ranges between K = 0.4 and 
K = 1.0. By contrast, the anharmonicity strongly modi- 
fies the TVBS-I properties. It is responsible for a strong 
redshift of the band as well as for a decrease of the band- 
width. The TVBS-I binding energy reaches 15.3 cm -1 
whereas the bandwidth is reduced to 19.0 cm -1 . 

Figs. 2 displays the TVBS wave function ^ K\{m) cor- 
responding to the situation described in Figs. 1 . Figures 
2a and 2b show the zero wave vector TVBS-I wave func- 
tion for A — and A — 8 cm -1 , respectively. The same 
curves are presented in figures 2c and 2d but for a wave 
vector equal to K = 0.69. The TVBS-I wave function 
is described by open squares whereas the TVBS-II wave 
function is represented by open circles. When K — 0, the 
spectrum supports a single bound state, namely TVBS-I 
(Figs. 1). When A = (Fig. 2a), the TVBS-I wave 
function is maximum when m = and slightly decreases 
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FIG. 1: Two-vibron energy spectrum for 
T = 5 K, xi = 25 pN, X 3 = 50 pN and 
for (a) A = and (b) A = 8 cm" 1 . 
The spectrum is centered onto the cor- 
rected frequency 2a>o and corresponds to 
the two-vibron dispersion curves drawn 
in half of the first Brillouin zone of the 
helix, i.e. < K < 7r. 
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FIG. 2: TVBS wave function ^ KX {m) for T = 5 K, X i = 25 
pN, X3 = 50 pN. (a) K = and A = cm" 1 , (b) K = and 
,4 = 8 cm" 1 , (c) K = 0.69 and A = cm _1 (d) K = 0.69 and 
A = 8 cm" 1 . The TVBS-I wave function is described by open 
squares whereas the TVBS-II wave function is represented by 
open circles. 



as the separating distance m between the two vibrons in- 
creases. It exhibits oscillations and takes significant val- 
ues for m = 1 and m = 3 whereas it almost vanishes when 
m = 2. When the anharmonicity is set to A = 8 cm -1 
(Fig. 2b), the maximum in m = is reinforced whereas 
the extension of the wave function has been reduced. In 
other words, such a state refers to two vibrons trapped 
around the same amide-I vibration. However, whatever 
the anharmonicity, the wave function shows a 3D charac- 
ter so that the bound state extends over the three spines 
of the helix. When K = 0.69, Figs. 1 clearly shows that 
the energy spectrum supports both TVBS-I and TVBS- 
II. The TVBS-I wave function appears strongly localized 
onto m — and exhibits basically the same features as 
when K = (Figs. 2c and 2d). By contrast, the TVBS- 



II wave function vanishes in m = 0. It is maximum in 
m = 3 and decreases as the separating distance between 
the two vibrons increases or decreases from m = 3 (Figs. 
2c and 2d). Such a state refers to two vibrons trapped 
around two nearest neighbor amide-I vibrations in a given 
spine. Nevertheless, the TVBS-II wave function shows a 
3D character so that the three spines of the helix are 
involved in the extension of the bound state. 

In Figs. 3, the two-vibron energy spectrum is illus- 
trated for T = 310 K and for the same set of parameters 
used in Figs. 1. Whatever the anharmonicity, the in- 
crease of the temperature strongly modifies the nature 
of the TVFS continuum. It appears symmetrically dis- 
tributed around 2co with a bandwidth which has been 
strongly reduced. It is equal to 33.4 cm -1 when A = 
(Fig. 3a) whereas it reaches 30.5 cm -1 when A = 8 
cm" 1 (Fig. 3b). As in Figs. 1, the TVBS-I lies be- 
low the continuum over the entire Brillouin zone. How- 
ever, its binding energy increases at biological tempera- 
ture whereas the corresponding bandwidth decreases. In 
the harmonic situation, i.e. A — (Fig. 3a), the TVBS-I 
binding energy is equal to 9.7 cm -1 and the bandwidth 
is equal to 8.4 cm" 1 . When the anharmonicity is fixed to 
A = 8 cm -1 (Fig. 3b), the binding energy reaches 25.9 
cm -1 whereas the bandwidth decreases to 3.6 cm -1 . As 
shown in Figs. 3, the spectrum displays now two wave 
vector windows for the occurrence of the TVBS-II. For 
A = (Fig. 3a), a first wave vector window appears in 
the range K — 0.7 and K = 1.4 and a second window 
takes place at the end of the first Brillouin zone, i.e. for 
K > 2.7. As previously, as when the anharmonicity is 
increased, the range of these two windows increases (Fig. 
3b) so that the TVBS-II band clearly tends to get out of 
the continuum. 

The corresponding TVBS wave functions are illus- 
trated in Figs. 4 for A = (Figs. 4a and 4c) and 
A = 8 cm" 1 (Fig. 4b and 4d). When K = (Figs. 4a 
and 4b), the spectrum supports the TVBS-I band, only. 
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FIG. 4: TVBS wave function t> K \(m) for T = 310 K, X i = 25 
pN, xs = 50 pN. (a) K = and A = cm -1 , (b) K = and 
,4 = 8 cnT 1 , (c) K = tt/3 and A = cm _1 (d) K = tt/3 and 
A — 8 cm _1 .The TVBS-I wave function is described by open 
squares whereas the TVBS-II wave function is represented by 
open circles. 



Whatever the anharmonicity, the TVBS-I wave function 
is maximum in to = and takes significant values for 
m = 3, m = 6, m = 9 ... and so on. When the anhar- 
monicity increases, the maximum in m = is reinforced 
and the extension of the wave function is reduced. In 
other words, at biological temperature, the two vibrons 
are confined in a given spine and are trapped around 
the same amide-I mode. When K — tt/3 (Figs. 4c and 
4d), Figs. 3 clearly show that the spectrum supports 
both TVBS-I and TVBS-II. The TVBS-I wave function 
is strongly localized onto to = and almost vanishes for 
the other m values. By contrast, the TVBS-II wave func- 
tion is clearly localized onto to = 3 so that it refers to 
two vibrons strongly trapped onto two nearest neighbors 
amide-I modes belonging to the same spine. Note that 



the anharmonicity does not modify significantly the na- 
ture of the wave function of both TVBS-I and TVBS-II. 

Finally, the two-vibron energy spectrum is illustrated 
in Figs. 5 for T = 310 K, X i = 37.5 pN and X 3 = 75.0 pN 
(E B (0) = 20.0 cm" 1 ). In that case, the strong vibron- 
phonon coupling is responsible for a drastic decrease of 
the TVFS bandwidth. As in Figs. 3, the anharmonic- 
ity enhances the decrease of the bandwidth which varies 
from 13.3 cm -1 for A = (Fig. 5a) to 11.2 cm -1 for 
A = 8 cm -1 (Fig. 5b). As previously, the spectrum ex- 
hibits a TVBS-I band located below the continuum. The 
anharmonicity reduces the corresponding bandwidth and 
yields a strong rcdshift of the band. However, the main 
difference when compared with the two previous cases is 
twofold. First, the dressing effect is strong enough to in- 
duce the occurrence of the TVBS-II band over the entire 
Brillouin zone, even when A = (Fig. 5a). The an- 
harmonicity increases the TVBS-II binding energy and 
reduces its bandwidth. Then, a third bound state band, 
called TVBS-III, occurs below the TVFS continuum. Al- 
though it does not appear clearly on the figure, this band 
lies below the continuum over the entire Brillouin zone. 
In a marked contrast with the two previous bound states, 
the TVBS-III binding energy decreases as when the an- 
harmonicity increases since it varies from 0.08 cm -1 when 
A = (Fig. 5a) to 0.076 cm" 1 when A = 8 cm" 1 (Fig. 
5b). By contrast, the TVBS-III bandwidth decreases 
with A. It is equal to 5.59 cm^ 1 when A — (Fig. 5a) 
and reaches 4.94 cm -1 when A = 8 cm -1 (Fig. 5b). 

In the latter situation corresponding to a strong dress- 
ing effect, we have verified that for both the TVBS-I and 
the TVBS-II, the two vibrons are confined in a given 
spine. TVBS-I refers to their trapping onto the same 
amide-I mode (to = 0) whereas TVBS-II characterizes 
their trapping onto two nearest neighbor amide-I modes 
(to = 3) located in the same spine (the wave functions 
are not drawn). The behavior of the wave function of 
the third bound state is illustrated in Figs. 6 for A = 
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FIG. 5: Two-vibron energy spectrum for 
T = 310 K, xi = 37.5 pN, xs = 75 pN 
and for (a) A = and (b) A = 8 cm -1 . 
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FIG. 6: TVBS-III wave function f K \(m) for T = 310 K, 
Xi = 37.5 pN, X3 = 75 pN. (a) K = tt/4 and A = cm" 1 , (b) 
K = tt/4 and A = 8 cm -1 , (c) K = tt/3 and A = cm _1 (d) 
K = tt/3 and A = 8 cm" 1 . 



(Figs. 6a and 6c) and A = 8 cm 1 (Figs. 6b and 6d). 
When K = tt/4 (Figs. 6a and 6b), the TVBS-III wave 
function exhibits non vanishing values for m = 1, 2, 4, 5, 
... and so on, with a resulting amplitude which decreases 
as m increases. By contrast, when K = tt/3 (Figs. 6c 
and 6d), the wave function appears strongly localized 
onto m = 1, only. These features indicate that TVBS- 
III refers to vibrons trapped onto neighboring amide-I 
modes belonging to two different spines. More precisely, 
when the TVBS-III binding energy is rather weak, for in- 
stance when K = 7r/4, the bound state is extended and 
the two vibrons are trapped over a large distance but be- 
long to two different spines. By contrast, the TVBS-III 
tends to localized when its binding energy increases so 
that it refers to vibrons trapped onto two nearest neigh- 
bor amidc-I modes belonging to two different spines. 



IV. DISCUSSION 

To interpret and discuss the previous results, let us 
first focus our attention onto the influence of the nonlin- 
ear sources, i.e. the intramolecular anharmonicity and 
the strong vibron-phonon coupling. As shown in Figs. 
1, 3 and 5, the intramolecular anharmonicity is respon- 
sible for a decrease of the TVFS bandwidth, typically of 
about a few wave numbers. This feature originates in 
the modification of the hopping constants due to both 
the anharmonicity and the dressing effect. Indeed, a 
TVFS corresponds to two independent vibrons moving 
according to the effective hopping constants J e //(n) = 
Ji(n)exp[— (1 + Arj)S{n,T)\ which differ from the con- 
stants involved in the harmonic approximation. Ba- 
sically, the intramolecular anharmonicity reinforces the 
role played by the coupling constant S(n, T) which char- 
acterizes the dressing effect (S(n, T) — > (1 + 4r])S(n, T)). 
Therefore, it enhances the dressing effect and favors a 
decrease of the effective hopping constants so that an- 
harmonic vibrons are more sensitive to the dressing ef- 
fect than harmonic vibrons. Note that as when the tem- 
perature is increased, the strength of the dressing effect 
increases which results in a drastic decrease of the TVFS 
bandwidth. As shown in Figs. 1, 3 and 5, the TVBS-I 
are located below the TVFS continuum over the entire 
Brillouin zone whereas, for TVBS-II, two situations oc- 
cur depending on the strength of the vibron-phonon cou- 
pling. For a weak coupling, the band disappears inside 
the continuum when the wave vector lies outside specific 
windows whereas for a strong coupling, the band is lo- 
cated below the continuum over the entire Brillouin zone. 
As when increasing the anharmonicity, the TVBS-I band 
is redshifted and its bandwidth is strongly reduced. In 
the same way, the anharmonicity modifies the nature of 
the TVBS-II band. If the band is resonant with the con- 
tinuum, the anharmonicity tends to open the wave vec- 
tor windows and when the TVBS-II band is not resonant 
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with the continuum, the anharmonicity induces a redshift 
of the band as well as a decrease of its bandwidth. Note 
that these different features generalize the work detailed 
in Ref. fljil ] devoted to the two-vibron dynamics in one- 
dimensional model of helix. We thus can conclude that 
the occurrence of both TVBS-I and TVBS-II originates 
in the interplay between the intramolecular anharmonic- 
ity A and the site dependent small polaron binding en- 
ergy Eb{w). More precisely, in a 3D helix, it is straight- 
forward to show that the small polaron binding energy 
takes significant values for n — and n — 3 and satisfies 
Eb(0) ~ 2Eb{3) in a perfect agreement with the results 
obtained in Ref. 0|. Therefore, A and Eb(0) are the key 
parameters which favors the occurrence of the TVBS-I 
through the terms Ab^b^bnbn in Eq.© which character- 
ize an attractive interaction between two vibrons located 
onto the same amide-I mode. By contrast, the occur- 
rence of TVBS-II is essentially controlled by the parame- 
ter Eb(3) which, as shown in Eq.JSJ), yields an attractive 
interaction between two vibrons located onto two near- 
est neighbor amide-I modes belonging to the same spine. 
To obtain a TVBS-II band below the TVFS continuum 
over the entire Brillouin zone, Eb(3) must exceed a crit- 
ical value typically of about 2J e //(3). As a consequence, 
TVBS-I refers to the trapping of the two vibrons onto 
the same amide-I mode whereas TVBS-II characterizes 
two vibrons trapped onto two nearest neighbor amide-I 
modes located in the same spine. 

However, the main results of the present study con- 
cerns the influence of the temperature on the nature of 
the two-vibron states in a real helix. 

As detailed in paper I , the single- vibron dynamics 
in a 3D a-helix is characterized by two kinds of hop- 
ping processes. The first kind refers to vibrational tran- 
sition between different spines of hydrogen-bonded pep- 
tide units whereas the second kind involves vibron hops 
along the same spine. In that context, the nature of 
the vibron dynamics essentially results from the compe- 
tition between these two kinds of processes. This previ- 
ous study has revealed that several parameters such as 
the temperature, the small polaron binding energy and 
the helix backbone conformation, allow for a transition 
between two regimes. At low temperature the small po- 
laron behaves like an undressed vibron so that the 3D 
nature of the helix plays a crucial role and the vibron is 
delocalized between the different spines. Within this un- 
dressed limit, vibron transitions between nearest neigh- 
bor and third nearest neighbor sites represent the dom- 
inant hopping mechanism since J(l) = —12.4 cm -1 and 
J(3) = 7.8 cm" 1 whereas J(2) = 3.9 cm -1 . By contrast, 
at biological temperature or strong small polaron bind- 
ing energy, the dressing effect strongly reduces the vi- 
brational exchanges between different spines so that the 
polaron tends to propagate along a single spine, only. 
Indeed, the dressing for vibrons hops between different 
spines is more efficient than the dressing which affects 
vibron transitions along the same spine so that J e //(3) 
becomes the dominant hopping constant. 



The singular behavior of the single vibron dynamics 
is responsible for the occurrence of a strong tempera- 
ture dependence on the two-vibron states. To illustrate 
this feature, let us first consider the temperature depen- 
dence of the TVFS continuum. Indeed, Figs. 1 clearly 
show that at low temperature the continuum exhibits a 
large bandwidth characterized by a complicated struc- 
ture whereas at biological temperature (Figs. 3 and 5), 
it appears symmetrically distributed around the central 
frequency 2lSq. In that context, it is straightforward to 
show that the continuum is perfectly reproduced from the 
knowledge of the single vibron energy spectrum and that 
the TVFS bandwidth is twice the single vibron band- 
width. Therefore, since single vibron states are delocal- 
ized over the three spines at low temperature, the dif- 
ferent hopping constants contribute to the bandwidth. 
For instance, when T = 5 K, \i — 25 pN, X3 — 50 pN 
and A — 8 cm -1 , the single vibron bandwidth is about 
50 cm" 1 and yields a TVFS bandwidth equal to 100.0 
cnr 1 in a perfect agreement with the results displayed 
in Figs. 1. Note that in that case the bandwidth mainly 
depends on both J e //(1) and J e //(3). By contrast, at bi- 
ological temperature, the continuum displayed in Figs. 3 
and 5 mainly refers to two independent vibrons confined 
in a given spine. When long range hopping constants 
are disregarded, the single vibron bandwidth is equal to 
4J e// (3). Therefore, when T = 310 K, X i = 25 pN and 
Xz = 50 pN, this value leads to a TVFS bandwidth equal 
to 29.12 cm" 1 when A = 8 cm -1 , in a perfect agree- 
ment with the value characterizing the real continuum 
(see Figs. 3). Note that the shape of the continuum 
results form the unfolding of the two-vibron dispersion 
curves of the ID model since the Brillouin zone of the 
helix is three time greater than the Brillouin zone of the 
corresponding ID model. Moreover, long range effects in 
the hop ping constants, which have been disregarded in 
Ref.|19j, lead to non vanishing values of the continuum 
for both K = ir/3 and K = n. 

The influence of the temperature on the two-vibron 
bound states is twofold. First, it strongly modifies the 
nature of both TVBS-I and TVBS-II. Then, it allows for 
the occurrence of a third bound state, i.e. TVBS-III, 
which has not been obtained within our previous one- 
dimensional model |19|. 

At low temperature (see Figs. 2), non vanishing inter- 
spine hopping constants yield a three dimensional nature 
of both TVBS-I and TVBS-II which the wave functions 
extend over the three spines of the helix. Although the 
two vibrons are trapped close to each other, there are not 
confined in the same spine so that the pair is allowed to 
explore the entire 3D nature of the helix. By contrast, 
at biological temperature, the dressing effect strongly re- 
duces the vibrational exchanges between different spines 
so that only intraspinc hopping constants remain. The 
pairs are thus confined in a given spine and exhibits the 
same features as the bound states described within the 
one-dimensional model (see Figs. 4). In a marked con- 
trast with the previous situation, the occurrence of the 
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third bound states is strictly due to the 3D nature of 
the helix. More precisely, as shown in Sec. II, this 3D 
nature manifests itself by the introduction of a general- 
ized site dependent small polaron binding energy Eb(ti). 
Although Eb(ji) takes significant values for n — and 
n = 3, it does not strictly vanish for n = 1. Note that 
J5s(n) ~ for n = 2 and n > 3. As a consequence, 
Eq.JSJ clearly shows that this non vanishing value con- 
tributes to an attractive interaction between two vibrons 
located onto two different spines. However, the formation 
of TVBS-III results from the competition between the 
interspine hopping constants, which tends to break the 
pair, and the previous attractive interaction which allows 
for the trapping of the pair. At low temperature, the in- 
terspine hopping constants are strong enough to prevent 
the occurrence of the third bound states. By contrast, 
at biological temperature, the dressing effect strongly re- 
duces the vibrational exchanges between different spines 
so that the attractive interaction mediated by Eg(l) is 
sufficient to create the trapping of the two vibrons onto 
two different spines. Note that strictly speaking the key 
parameter for the occurrence of the TVBS-III is B(l), 
which, in addition to its strong dependence on _Eb(1), 
exhibits a weak contribution B(l) proportional to the 
intramolecular anharmonicity (see Eq.JJJJO- Due to the 
sign of the hopping constant J(l), B(l) is negative and 
tends to produce a repulsive interaction between the two 
vibrons. As a consequence, the binding energy of TVBS- 
III slightly decreases with the anharmonicity as shown in 
Figs. 5. 

At this step, let us mention that this theoretical 
work clearly suggests the occurrence of two-vibron bound 
states for the amide-I vibrations in a-helices. In a recent 
experiment, P. Hamm and co-workers [2lL |22| have per- 
formed the infrared pump-probe spectroscopy of the N-H 
stretch vibrations in a stable a-helix, the poly-7-benzyl- 
L-glutamate (PBLG). They observed two excited-state 
absorption bands which were interpreted in terms of the 
occurrence of two kind of two-vibron bound states. In 
that context, the knowledge of both the single- vibron 
and the two-vibron dynamics allows us to simulate this 
pump-probe spectrum but for amide-I vibrations. Dur- 
ing a pump-probe experiment, the helix interacts with 
two infrared short pulses, i.e. the pump and the probe 
pulse HI El HJ |H IH H3. The pump pulse excites 
the helix from its vibrational ground state to the first ex- 
cited state. In other words, the pump creates a zero wave 
vector single vibron state. By contrast, the probe pulse 
measures the change of absorption as a result of that exci- 
tation. As a consequence, the pump-probe spectrum ex- 
hibits three contributions. The first contribution, called 
the Bleach, corresponds to the loss of absorption due to 
the depletion of the ground state whereas the second con- 
tribution characterizes the stimulated emission from the 
first excited state. These two processes lead to a negative 
peak in the pump-probe signal. The third contribution, 
which accounts for the creation of zero wave vector two- 
vibron states from the excited single-vibron state, yields 
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FIG. 7: Pump-probe spectrum of the amide-I band for (a) 
T = 310 K, xi = 25 pN, xa = 50 pN, A = 8 cm" 1 and (b) 
T = 310 K, xi = 37.5 pN, X s = 75 pN, A = 8 cm" 1 . 



a positive peak in the spectrum. From the knowledge 
of both the single-vibron 0] and two-vibron states, it 
is straightforward to calculate the matrix element of the 
dipole involved in the three transitions characterizing the 
spectrum |33j . To proceed, both the pump and the probe 
pulses are assumed to be polarized along the axis of the 
helix and an artificial width of 1 cm" 1 has been intro- 
duced to plot the spectrum. 

Fig. 7a shows the pump-probe bare spectrum at T = 
310 K and with the parameter values as in Fig. 3b. The 
spectrum, which is centered onto 2cJo, exhibits a single 
negative peak located at -7.2 cm -1 and a single positive 
peak located at -34.0 cm -1 . In agreement with the result 
displayed in Fig. 3b, the positive peak is the signature of 
the zero wave vector bound state TVBS-I supported by 
the helix. Note that the frequency difference between the 
two peaks is about the binding energy of TVBS-I. For 
a stronger vibron-phonon coupling, i.e. for parameter 
values as in Fig. 5b, the pump-probe spectrum shows 
three peaks. A negative peak is located at -2.7 cm" 1 
whereas two positive peaks occur at -18.6 cm" 1 and - 
62.6 cm -1 , respectively. As shown in Fig. 5b, the two 
positive peaks characterize the two bound states TVBS- 
I and TVBS-II. The intensity of the peak connected to 
TVBS-II is about two times greater than the intensity 
of the peak associated to TVBS-I. Note that the binding 
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energy of TVBS-III is too weak so that the corresponding 
peak is embedded in the bleach. 

As illustrated in Figs. 7, a pump-probe experiment 
provides a direct evaluation of the strength of the differ- 
ent nonlinear sources favoring the occurrence of bound 
states through the measurement of both the number of 
positive peaks and the values of the corresponding bind- 
ing energy. Nevertheless, let us mention that relaxation 
mechanism have been neglected in the present analysis 
so that the width of the different peaks has been disre- 
garded. For amide-I vibration, we expect a competition 
between the binding energy of the bound states and the 
width of the corresponding peaks which will result in a 
difficulty to separate the different peaks. Note that it 
was not the case for N-H modes in PGLB because al- 
though the width of the different peaks was about 50 
cm" 1 , the binding energies of both TVBS-I and TVBS-II 
were about 275 and 120 cm -1 , respectively. These large 
values originate in the strong anharmonicity of each N-H 
stretching mode (A — 60 cm -1 ) as well as in the strong 
vibron-phonon coupling (i£g(0) ~ 84 cm" 1 ). 

To conclude this discussion, let us address the funda- 
mental question of the breather-like behavior of the two- 
vibron bound states. Indeed, the formation of discrete 
breathers in nonlinear classic lattices has been the sub- 
ject of intense theoretical research during the last decade 
(for a recent review, see for instance Refs. |l|, |2J, y, |8| ) . 
Discrete breathers, which correspond to highly localized 
vibrational excitations in anharmonic lattices, do not re- 
quire integrability for their existence and stability. They 
are not restricted to one-dimensional lattices and it has 
been suggested that they should correspond to quite gen- 
eral and robust time periodic solutions. Recently, the 
quantum nature of discrete breathers has been investi- 
gated and it has been suggested that bound states involv- 
ing several quanta are good candidates for the quantum 
counterpart of discrete breathers 0, 0, IE El IE IS HE • 

To illustrate the concept of quantum breather in a- 
helices, we restrict our attention to TVBS-I which are 
located below the TVFS continuum over the entire first 
Brillouin zone. As shown previously, these states cor- 
respond to the trapping of two vibrons onto the same 
amide-I mode and described two quanta which behave as 
a single particle. Note that, at this step, the breather-like 
behavior manifests itself in the localization of the sepa- 
rating distance between the two vibrons which form the 
pair. To characterize the motion of the pair, we use the 
procedure introduced by Proville [3(j which consists in 
defining the time dependent Wannier wave function as 

v A K 

When A is restricted to the TVBS-I band, <fi n (t) can 
be viewed as the wave function associated to the two 
trapped vibrons viewed as a single particle. To charac- 
terize the motion of that particle, we introduce a coher- 
ent lifetime, proportional to the invert of the TVBS-I 
bandwidth, which measures the time spent by the pair 
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FIG. 8: Evolution of the probability P„(t) to observe the 
vibron pair on the site n at time t. The pair, which has been 
created on the site n = at time t = 0, is restricted to the 
TVBS-I band (see the text). (a) T = 5 K, X i = 25 pN, X 3 = 50 
pN, A = 8 cm" 1 , (b) T = 310 K, X i = 25 pN, X3 = 50 pN, 
A = 8 cm" 1 , (c) T = 310 K, X i = 37.5 pN, X3 = 75 pN, 
A = 8 cm" 1 . 



on a given site. A zero bandwidth produces an infinite 
coherent lifetime for the pair which exhibits a breather- 
like behavior, i.e. the two vibrons are trapped around 
the site where the pair has been created. By contrast, 
a finite bandwidth leads to the derealization of the pair 
which moves along the lattice. 

By assuming that the pair is created onto the site n — 
at the initial time t = 0, Figs. 8 show the behavior of 
the probability P n (t) =| <j> n (t) | 2 to observe the pair onto 
the site n at time t. The three situations correspond- 
ing to the parameters used in Fig. lb, 3b and 5b are 
illustrated in Fig. 8a, 8b and 8c, respectively. At low 
temperature, Fig. 8a clearly shows the 3D nature of the 
TVBS-I since the pair tends to rapidly delocalize over the 
three spines of the helix. Indeed, after a coherent lifetime 
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of about 1.5 ps, 80 % of the initial population have left 
the excited site n — and have been tranfered to sites 
belonging to the other spines. Note that the population 
of the excited site shows damped oscillations due to the 
strong anistropy which discrimites between inter and in- 
tra spine vibrational exchanges. When the temperature 
reaches T — 310 K, the confinement of the pair takes 
place as illustrated in Fig. 8b so that the pair propagates 
along the spine where it has been created. Nevertheless, 
the finite value of the TVBS-I bandwidth prevents the 
localization of the pair which leaves the excited site af- 
ter a coherent lifetime of about 7.8 ps. However, for a 
stronger vibron-phonon coupling, the TVBS-I bandwidth 
is drastically reduced so that the breather-like behavior 
of the pair manifests itself by an enhancement of its local 
character (Fig. 8c). As previously, the pair is confined 
in the spine where it has been created but its coherent 
lifetime has been drastically increases and is about 86 ps. 

V. CONCLUSION 

In the present paper, the two-vibron dynamics associ- 
ated to amide-I modes in an a-helix has been described. 
The 3D nature of the helix was modeled by three spines 
of hydrogen-bonded peptide units linked via covalent 
bonds. To remove the intramolecular anharmonicity of 
each amide-I mode and to rcnormalizc the strong cou- 
pling between the vibrons and the phonons associated to 
the external motion of the residues, two unitary transfor- 
mations have been applied to finally obtain the dressed 
anharmonic vibrons point of view. In a perfect agree- 
ment with recent calculations performed within a one- 
dimensional model of helix, it has been shown that the 
two-vibron energy spectrum supports both a TVFS con- 
tinuum and two kinds of bound states, called TVBS-I 
and TVBS-II, connected to the trapping of two vibrons 
onto the same amide-I mode and onto two nearest neigh- 
bor amide-I modes belonging to the same spine, respec- 
tively. In addition, it has been shown that the interplay 



between the temperature and the 3D nature of the he- 
lix strongly affects the characteristics of the two-vibron 
states. This feature originates in the singular behavior 
of the single vibron dynamics which displays a transition 
between two regimes. At low temperature, the single 
polaron behaves as an undressed vibron delocalized be- 
tween the different spines whereas at biological tempera- 
ture, the dressing effect strongly reduces the vibrational 
exchanges between different spines so that the polaron is 
confined in a single spine. The consequence is twofold. 
First, the nature of both TVBS-I and TVBS-II is strongly 
modified. Then, a third bound state called TVBS-III oc- 
curs. At low temperature, non vanishing interspine hop- 
ping constants yield a three dimensional nature of both 
TVBS-I and TVBS-II which the wave functions extend 
over the three spines of the helix whereas, at biologi- 
cal temperature, the pairs are confined in a given spine 
and exhibits the same features as the bound states de- 
scribed within a ID model. The occurrence of the TVBS- 
III results from the competition between the interspine 
hopping constants and the site dependent small polaron 
binding energy which favors an attractive interaction be- 
tween vibrons located onto different spines. If at low 
temperature the hopping constants are strong enough to 
prevent the occurrence of TVBS-III, it is no longer the 
case at biological. As a result the attractive interaction 
becomes the dominant mechanism leading to the trap- 
ping of the two vibrons onto two different spines. Fi- 
nally, the breather-like nature of the TVBS-I has been 
analyzed by studying the time evolution of a vibron pair 
initially created on a given site. Due to the 3D nature of 
TVBS-I, the pair rapidly delocalizes over the three spines 
at low temperature which results in a coherent lifetime 
of about a few picoseconds. By contrast, at biological 
temperature, the pair is confined in the spine where it 
has been created. Its coherent lifetime increases with the 
strength of the vibron-phonon coupling and it has been 
shown that a strong coupling can localize the pair during 
a time of about 86 ps. 
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